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We derive general conditions for 100% frequency conversion in any doubly resonant nonlinear cav- 
ity, for both second- and third- harmonic generation via x^'^^ ^ind x^^^ nonlinearities. We find that 
conversion efficiency is optimized for a certain "critical" power depending on the cavity parameters, 
and assuming reasonable parameters we predict 100% conversion using milliwatts of power or less. 
These results follow from a semi-analytical coupled-mode theory framework which is generalized 
from previous work to include both x^'^^ X^^^ media as well as inhomogeneous (fully vectorial) 
cavities, analyzed in the high-efficiency limit where down- conversion processes lead to a maximum 
efficiency at the critical power, and which is verified by direct finite-difference time-domain (FDTD) 
simulations of the nonlinear Maxwell equations. Explicit formulas for the nonlinear coupling coeffi- 
cients are derived in terms of the linear cavity eigenmodes, which can be used to design and evaluate 
cavities in arbitrary geometries. 



PACS numbers: 

INTRODUCTION 

In this paper, we consider second- and third-harmonic 
generation in doubly resonant cavities. We generalize 
previous experimental and theoretical work on this sub- 
ject, which had focused only on x^^^ nonlinearities in 
large Fabry-Perot etalons or two-dimensional ring res- 
onators where a scalar approximation applied [U [2l |3l 
H El El O E], to incorporate both x^^^ and x^^^ non- 
linearities and handle the inhomogenous fully vectorial 
case. We then develop several results from this gen- 
eralization: whereas it is well known that 100% har- 
monic conversion is possible, at least in x^^^ media, in 
this work we further explore and identify the conditions 
under which this can be achieved. First, we demon- 
strate the existence of a critical input power at which 
harmonic generation is maximized, in contrast to previ- 
ous work that focused largely on the low-power limit in 
which generation efficiency increased monotonically with 
input power [3^. Second, while it is well known that har- 
monic conversion can be achieved at arbitarily low pow- 
ers given sufficiently long cavity lifetimes, this implies a 
narrow bandwidth — we show that, by combining mod- 
erate lifetimes (0.1% bandwidth) with tight spatial con- 
finement, 100% second- and third-harmonic conversion 
can theoretically be achieved with sub-milliwatt power 
levels. Such low-power conversion could find applica- 
tions such as high-frequency sources [5], ultracompact- 
coherent optical sources [TOl [Til HI], imaging [1^, and 
spectroscopy |14 . 

Nonlinear frequency conversion has been commonly re- 
alized in the context of waveguides [151 [ISl [III US? oi" 
even for free propagation in the nonlinear materials, in 
which light at one frequency co-propagates with the gen- 
erated light at the harmonic frequency |20l [211 ^ • 
A phase-matching condition between the two frequen- 
cies must be satisfied in this case in order to obtain ef- 



ficient conversion [5, 9 . Moreover, as the input power 
is increased, the frequency conversion eventually satu- 
rates due to competition between up and down conver- 
sion. Frequency conversion in a doubly resonant cavity 
has three fundamental differences from this familiar case 
of propagating modes. First, light in a cavity can be 
much more intense for the same input power, because 
of the spatial (modal volume V) and temporal (lifetime 
Q) confinement. We show that this enhances second- 
harmonic (x^^^) conversion by a factor of /V and en- 
hances third-harmonic (x^^^) conversion by a factor of 
/V . Second, there is no phase-matching condition per 
se for 100% conversion; the only absolute requirement 
is that the cavity support two modes of the requisite 
frequencies. However, there is a constant factor in the 
power that is determined by an overlap integral between 
the mode field patterns; in the limit of a very large cavity, 
this overlap integral recovers the phase-matching condi- 
tion for x^^^ processes. Third, the frequency conversion 
no longer saturates — instead, it peaks (at 100%, with 
proper design) for a certain critical input power satisfy- 
ing a resonant condition, and goes to zero if the power is 
either too small or too large. 

Second-harmonic generation in cavities with a single 
resonant mode at the pump frequency [lOl HH [23l [24l 

25, 261 ^ ^ ^ [Ml [311 [32 133 [Ml [33 [Ml [33 [sg 

or the harmonic frequency [39] requires much higher 
power than a doubly resonant cavity, approaching one 
Watt [3 |38] and/or requiring amplification within the 
cavity. (A closely related case is that of sum-frequency 
generation in a cavity resonant at the two frequencies 
being summed [TO.) Second-harmonic generation in a 
doubly resonant cavity, with a resonance at both the 
pump and harmonic frequencies, has most commonly 
been analyzed in the low-efficiency limit where nonlin- 
ear down- conversion can be neglected [H [3 El [3 El [9] , 
but down-conversion has also been included by some au- 
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thors m m [3]. Here, we show that not only is down- 
conversion impossible to neglect at high conversion effi- 
ciencies (and is, in fact, necessary to conserve energy), 
but also that it leads to a critical power where harmonic 
conversion is maximized. This critical power was demon- 
strated numerically by Ref. |4T] in a sub-optimal geom- 
etry where 100% efficiency is impossible, but does not 
seem to have been clearly explained theoretically; the 
phenomenon (for x^^^) was also implicit in the equations 
of Ref. 3 but was not identified, probably because it oc- 
curred just beyond the range of power considered in that 
work. 
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FIG. 1: Top: Schematic diagram of waveguide-cavity system. 
Input light from a waveguide (left) at one frequency (ampli- 
tude si+) is coupled to a cavity mode (amplitude ai), con- 
verted to a cavity mode at another frequency (amplitude 02) 
by a nonlinear process, and radiated back into the waveguide 
(amplitude S2-). Reflections at the first frequency (si-) may 
also occur. Bottom: Id example, formed by quarter- wave de- 
fect in a quarter-wave dielectric stack. Dielectric material is 
yellow, and electric field Ez of third-harmonic mode is shown 
as blue/ white/red for positive/zero/negative amplitude. 

Previous work on third-harmonic generation in cavi- 
ties considered only singly resonant cavities; moreover, 
past work focused on the case of x^^^ materials where Scj 
is generated by cascading two nonlinear processes (har- 
monic generation and frequency summing) [141 142]. Here, 
we examine third-harmonic generation using x^^-^ mate- 
rials so that only a single resonant process need be de- 
signed and a different set of materials becomes available, 
(x^^-* third-harmonic generation in a bulk periodic struc- 
ture, with no cavity, was considered in Ref. |43l) In a 
^(3) medium, there are also self/cross-phase modulation 
phenomena (nonlinear frequency shifts) that, unchecked, 
will prevent 100% conversion by making the frequency 
ratio 7^ 3. To address this mismatch, we describe how 
one can use two materials with opposite-sign x^^^ to can- 
cel the frequency-shifting effect; it may also be possible 
to pre-shift the cavity resonant frequency to correct for 
the nonlinear shift. On the other hand, a x^^^ medium 
has no self-phase modulation, and so in this case it is suf- 
ficient to increase the input power until 100% frequency 
conversion is reached. (An "effective" self-phase mod- 
ulation occurs in x*^^^ media due to cascaded up- and 
down- conversion processes [44 , but these processes are 
fully taken into account by our model. We also consider 
media with simultaneous x^^^ and x^^^ nonlinear it ies, and 
show that the latter can be made negligible.) If the criti- 
cal field were too intense, then material breakdown might 



also be an obstacle, but we show that it is sufficient to 
use modes with a large lifetime Q and small volume V 
so that a slow conversion due to a weak nonlinear effect 
has enough time to occur. 

In particular, we consider the general situation de- 
picted schematically in Fig.[l] a two- mode nonlinear cav- 
ity coupled to an input/output channel. For example, a 
one-dimensional realization of this is shown in Fig. [T] a 
Fabry-Perot cavity between two quarter- wave stacks [45] , 
where the stack has fewer layers on one side so that light 
can enter /escape. For a nonlinear effect, we consider 
specifically a x^^^ nonlinear ity, corresponding essentially 
to a shift in the refractive index proportional to the non- 
linear susceptibility x^^^ multiplied by electric field E to 
the (£ — 1)^^ power. Most commonly, one would have 
either a x^^^ (Pockels) or x^^^ (Kerr) effect. Such a non- 
linearity results in harmonic generation ^ : light with 
frequency uo is coupled to light with frequency lij. There- 
fore, if we design the cavity so that it supports two modes, 
one at uo and one at ^uo^ then input power at uo can be 
converted, at least partially, to output power at iuo. 

In the following, we derive a semi-analytical descrip- 
tion of harmonic generation using the framework of 
coupled-mode theory [Iliai3[5l[S[T0l[ll|3ll[3a|47], 
and then check it via direct numerical simulation of the 
nonlinear Maxwell equations [71 |48l [49] . For maximum 
generality, we derive the coupled-mode equations using 
two complementary approaches. First, we use "tempo- 
ral" coupled-mode theory [50l [51] , in which the general 
form of the equations is determined only from principles 
such as conservation of energy and reciprocity, indepen- 
dent of the specific physical problem (for example, elec- 
tromagnetic or acoustic waves). Second, we apply per- 
turbation theory directly to Maxwell's equations in order 
to obtain the same equations but with specific formulas 
for the coupling coefficients in terms of the linear eigen- 
modes. Unlike most previous treatments of this prob- 
lem [TJ [2j [47] , we do not make a one-dimensional or scalar 
approximation for the electromagnetic fields (invalid for 
wavelength-scale cavities), and we consider both x^^^ and 
^(3) media. (The optimization of these coupling coeffi- 
cients is then the generalization of the phase-matching 
criteria used in one-dimensional geometries [5].) 

TEMPORAL COUPLED-MODE THEORY 

We derive coupled-mode equations describing the in- 
teraction of light in a multi-mode cavity filled with non- 
linear material and coupled to input/output ports, from 
which light can couple in (s+) and out (5_) of the cav- 
ity. A schematic illustration of the system is shown in 
Fig. [1] Specifically, we follow the formalism described 
in Ref. 50, adapted to handle nonlinear ly coupled modes 
with frequencies uok- Although similar equations for the 
case of x^^^ media were derived in the past [J, they do 
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not seem to have been derived for x^^^ harmonic gener- 
ation in cavities. Moreover, a derivation via the tempo- 
ral coupled- mode formalism of Ref. |50] is arguably more 
general than earlier developments based on a particular 
scalar nonlinear wave equation, because this formalism 
(for a given-order nonlinearity) depends only on gen- 
eral considerations such as weak coupling and energy 
conservation (the resulting equations hold for vector or 
scalar waves in electromagnetism, acoustics, or any other 
weakly-coupled problem with a few simple properties). In 
the next section, we will then specialize the equations to 
electromagnetism by deriving explicit equations for the 
coupling coefficients from Maxwell's equations. 

We let a/e denote the time-dependent complex ampli- 
tude of the kth mode, normalized so that |a/cp is the 
electromagnetic energy stored in this mode. We let s± 
denote the time-dependent amplitude of the incoming 
(+) or outgoing (— ) wave, normalized so that is 
the power. (More precisely, s±{t) is normalized so that 
its Fourier transform |s±(a;)p is the power at co. Later, 
we will let Sk± denote the input /output power at cok-) 
[In Id, the units of |a/cp and \s±\'^ are those of energy 
and power per unit area, respectively. More generally, 
in d dimensions, the units of |a/cp and \s±\'^ are those of 
energy and power per length^" By itself, a linear cav- 
ity mode decaying with a lifetime Tk would be described 
by dak/dt = {iuok — l/rk)ak. [Technically, such a decay- 
ing mode is not a true eigenmode, but is rather a "leaky 
mode" [52j, corresponding to a "quasi-bound state" in 
the Breit-Wigner scattering theory [53j.] The decay rate 
I/t/c can be decomposed into l/r^ = l/re^k-^^/^s,k where 
l/Te,/c is the "external" loss rate (absorption etc.) and 
I/ts^/c is the decay rate into s_. When the weak coupling 
{uJkTk ^ 1) to s± is included, energy conservation and 
similar fundamental constraints lead to equations of the 
form |54.: 



dak f. l\ I 2 

dt \ TkJ V Ts^k 
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(2) 



This can be generalized to incorporate multiple in- 
put/output ports, direct coupling between the ports, and 
so on [51]. The only unknown parameters in this model 
are then the frequencies uo^ and the decay rates 1/r/e, 
which can be determined by any numerical method to 
solve for the cavity modes (e.g. FDTD, below). Instead 
of T/c, one commonly uses the quality factor Qk = oOkTk/2. 

Nonlinearity modifies this picture with two new 
amplitude-dependent effects: a shift in the frequency 
(and decay rate) of the cavity, and a coupling of one 
cavity mode to another. We neglect nonlinear effects on 
the input/output ports, under the assumption that in- 
tense fields are only present in the cavity (due to spatial 
and temporal confinement). We will also make two stan- 
dard assumptions of nonlinear systems. First, that the 
nonlinearities are weak, in the sense that we can neglect 
terms of order {x^^^Y higher; this is true in practice 
because nonlinear index shifts are always under 1% lest 
material breakdown occur. Second, we make the rotat- 
ing wave approximation: since the coupling is weak, we 
only include terms for ak that have frequency near uok- 
In particular, we suppose that uOk ^ kuoi^ so that uOk is 
the /cth harmonic. The result is that, for a given order 
nonlinearity, there are only a few possible new terms that 
can appear in the coupled-mode equations. In particular, 
for a x^^^ nonlinearity with two modes uji and its second 
harmonic 002^ the coupled- mode equations must take the 
form: 



dai 
~dt 

da 2 
~dt 



l\ 2 

iui I — iujipiala2 + w 5+ (3) 

iU2 - — ] a2 - iu2P2Ci\ + 4/^-5+ (4) 

T2 J V Ts^2 



Similarly, for a x^^^ nonlinearity with two modes uJi and 
its third harmonic ojs, the coupled- mode equations must 
take the form: 



dai 
~dtt 

das 
~dtt 



- ail \ai\^ - Q^i3 \as\^ ) - ai 
[l - ^33 \a2,\^ - ^31 |ai|^^ - 



iui(3i{a\) as^il s+ 




as - icJsPsa^ + w 5+ 

^3/ V Ts^S 




(5) 
(6) 



In equations [5I-I6) one sees two kinds of terms. The first are frequency -shifting terms, with coefficients a^j, depen- 
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dent on one of the field amplitudes. For x*^^\ this effect is 
known as self-phase and cross-phase modulation, which 
is absent for x*^^^ (under the first-order rotating- wave ap- 
proximation). The second kind of term transfers energy 
between the modes, with coupling coefficients /3i^ corre- 
sponding to four- wave mixing for x^^^- Furthermore, we 
can constrain the coupling terms Pi by energy conser- 
vation: -^{\aif + |<^2|^) = 0. For x*^^\ the constraint 
that follows is: coipi = 002^21 X^^\ the constraint 
is uJiPi = uJsPs' (This constraint holds even in cavities 
with external loss as discussed in Sec. : energy is still 
conserved in the sense that the input power must equal 
the output power plus the loss power, and so the har- 
monic conversion term must lead to an equal energy loss 
and gain at uji and uj2,3, respectively.) 

The general process for construction of these coupled- 
mode equations is as follows. The underlying nonlinear- 
ity must depend on the physical, real part of the fields, 
corresponding to {ak + <^fc)/2. It then follows that the 
^(^) term will have i powers of this real part, giving var- 
ious product terms like aja2 (for x*^^^) a^aiai (for 
X^^^). Most of these terms, however, can be eliminated 
by the rotating-wave approximation. In particular, we 
assume that each ak term is proportional to e^*'^ multi- 
plied by a slowly varying envelope, and we discard any 
product term whose total frequency differs from kuj for 
the dak/dt equation. Thus, a term like alasas would be 
proportional to e^*'^, and would only appear in a da^/dt 
equation. (We focus on the simpler case of doubly reso- 
nant cavities in this paper.) 

At this point, the equations are already useful in or- 
der to reason about what types of qualitative behaviors 
are possible in general. In fact, they are not even spe- 
cific to electromagnetism and would also apply to other 
situations such as acoustic resonators. However, in or- 
der to make quantitative predictions, one needs to know 
the nonlinear coefficients a^j and Pi (as well as the linear 
frequencies and decay rates). The evaluation of these co- 
efficients requires a more detailed analysis of Maxwell's 
equations as described below. 

PERTURBATION THEORY AND COUPLING 
COEFFICIENTS 

In this section, we derive explicit formulas for the non- 
linear coupling coefficients in the coupled-mode theory of 
the previous section, applied to the case of electromag- 
netism. Unlike previous work, our expressions apply to 



the fully vectorial equations, valid for high index-contrast 
materials, and we derive the x*^^^ case as well as x*^^^- Our 
derivation is closely related to that of Ref. which only 
considered the frequency shifting (self-phase modulation) 
and not harmonic generation. 

When a dielectric structure is perturbed by a small Se, 
a well-known result of perturbation theory states that 
the corresponding change Sco in an eigenfrequency co is, 
to first order [45] : 

6ij _ l /d^xfe|Ep _ l /d^xE*-(5P 
~"2 J^3xe|E|' ~"2 /d3x£|E|' 

where E is the unperturbed electric field and JP = feE 
is the change in polarization density due to Se. In fact, 
Eq.[7|is general enough to be used with any ^P, including 
the polarization that arises from a nonlinear susceptibil- 
ity. In particular, we can use it to obtain the coupling 
coefficients of the GMT. 

To do so, we first compute the nonlinear first-order fre- 
quency perturbation due to the total field E present from 
all of the modes. Once the frequency perturbations 6uJk 
are known, we can re-introduce these into the coupled- 
mode theory by simply setting cok ^ ^^k in Eq. [l] 
By comparison with Eqs. |3]-[6| we can then identify the 
a and p coefficients. 

We consider first a x^^^ nonlinearity, with the nonlinear 
polarization SP given by SPi = Y^^.^eXiji^jEk, in a 
cavity with two modes Ei and E2. As before, we require 
that the modes oscillate with frequency uoi and 002 ^ 2cJi , 
respectively. Taking E = Re[Eie''^i^ + E2e^'^2t] ^nd using 
the rotating-wave approximation, we can separate the 
contribution of JP to each Jcj/c, to obtain the following 
frequency perturbations: 

5u:, ^ 1 / ^^x Y.^Jk ^xgl [El, {E2jEl, + ^r,^2fe) ] 
^1 4 Jd^^e\Eif 

(8) 

5u;2 _ 1 / Eijfc ^Xjfk E^-EijEik 
002 ~ ^ jd^^e\E2\^ 

Similarly, for a centro-symmetric x^^^ medium, ^P is 
given by 8V = £x^^^|EpE, with E = Re[Eie^^i^ + 
Ese*'^^^]. We obtain the following frequency perturba- 
tions: 
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5uJi 



f fcx^') (|Ei • Ell V 2 |Ei + 2(Ei • Et)(E3 • E|) 



/d3x£|Ei| 



2|Ei ■E3|' + 2|Ei -E^l' 



■3(EI-EI)(Ef -Es)) 



/d3x£|Ei| 



(10) 



Su)3 
UJ3 



/d3x£x(3) (IE3 • E3I' + 2 IE3 • E||' + 2(Ei • Et)(E3 • E|) 



/rf3x£|E3 



+ 



2 |Ei ■ E3I' + 2 |Ei ■ E^l' + (El • Ei)(Ei • E^)) 



/d3x£|E3| 



(11) 



There is a subtlety in the apphcation of perturbation the- 
ory to decaying modes, such as those of a cavity coupled 
to output ports. In this case, the modes are not truly 
eigenmodes, but are rather "leaky modes" [52 , and are 
not normalizable. Perturbative methods in this context 
are discussed in more detail by [511 152] , but for a tightly 
confined cavity mode it is sufficient to simply ignore the 
small radiating field far away from the cavity. The field 
in the cavity is very nearly that of a true eigenmode of 
an isolated cavity. 

As stated above, we can arrive at the coupling coef- 
ficients by setting uj^ ^ ^^k in Eq. [l] However, 
the frequency perturbations 5uOk are time-independent 
quantities, and we need to connect them to the time- 
dependent a/e amplitudes. Therefore, to re-introduce the 
time dependence, one can use the slowly varying enve- 
lope approximation: a slowly varying, time-dependent 
amplitude ak{t) is introduced into the unperturbed fields 
E/e E/ca/c(t). The eigenmode must be normalized so 
that |a/ep is the energy, as assumed for the coupled- mode 



theory. Thus, we divide each E/^ by y ^ / £:|E/cp. 

First, we consider the x^^^ medium. Carrying out the 
above substitutions in Eq. [l] and grouping terms propor- 
tional ak yields Eqs. l3||4] with a^j and Pi given by: 



(12) 



1/2 



/?2 



^ [/d3x£|Ei|'] [/d3x£|E2|' 

1 /d'x Y.iWlE2iEi,Eik 



1/2 



(13) 
(14) 



/d3x£|Ei|' /d3x£|E2|' 

A similar calculation yields the x^^^ coupled-mode 



equations with coefficients given by: 



l/d3x £x(3)|Ei-Ei|V|Ei-E* 



/rf3x £|Ei 



(15) 



"31 



1 /d3x £x(3) |Ei|2 IE3I' + |Ei • E3I' + |Ei ■ E^ 



/cZ3x£|Ei|^ /d3x£|E: 



ai3 = asi 

3 /d3x£x(3)(EI-EI)(EI-E3) 



/d3x £|Ei| 



3/2 



/rf3x£|E3 



1/2 



/?3 = 



1 /d3x£x(3)(Ei-Ei)(Ei-E|) 



/d3x £|Ei| 



3/2 



/rf3x£|E3 



1/2 



(16) 
(17) 
(18) 



Note that Eqs. p!2||T8| verify the conditions uoijSi = 
002^2 and uoijSi = 00^,131^ previously derived from con- 
servation of energy — for x^^\ this requires that one ap- 
ply the symmetries of the Xi/k tensor, which is invariant 
under permutations of ijk for a frequency-independent 
^(2) [45]^ Furthermore, we can relate the coefficients a 
and j3 to an effective modal volume V ^ similar to Ref. l55l 
In particular, the strongest possible nonlinear coupling 
will occur if the eigenfields are a constant in the nonlin- 
ear material and zero elsewhere. In this case, any integral 
over the fields will simply yield the geometric volume V 
of the nonlinear material. Thus, for the x^^^ effect we 
would obtain pi ~ x^^V^/^5 similarly, for the x^^^ effect 
we would obtain aij^/3i ~ x*^^V^^- This proportionality 



to 1/VV and 1/V carries over to more realistic field pro- 
files (and in fact could be used to define a modal volume 
for these effects). 
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NUMERICAL VALIDATION 

To check the predictions of the x^^^ coupled- mode 
equations, we performed an FDTD simulation of the 
one-dimensional waveguide-cavity system shown in Fig[l] 
whose analytical properties are uniquely suited to third- 
harmonic generation. (The FDTD method, including 
techniques to simulate nonlinear media, is described 
in Ref. ^56,) This geometry consists of a semi-infinite 
photonic-crystal structure made of alternating layers of 
dielectric (^i = 13 and £2 = 1) with period a and 
thicknesses given by the quarter-wave condition {di = 
V^2/(V^i + V^2) ^2 = a — di^ respectively). Such a 
quarter-wave stack possesses a periodic sequence of pho- 
tonic band gaps centered on frequencies coi = (y^^ + 
V^2)/^-\/^i^2 (in units of 27rc/a) for the lowest gap, and 
higher-order gaps centered on odd multiples of ui . More- 
over, a defect formed by doubling the thickness of a £1 
layer creates cavity modes at exactly the middle of ev- 
ery one of these gaps. Therefore, it automatically satis- 
fies the frequency-matching condition for third-harmonic 
generation. In fact, it is too good: there will also be 
"ninth harmonic" generation from cjs to cjg. This un- 
wanted process is removed, however, by the discretization 
error of the FDTD simulation, which introduces numer- 
ical dispersion that shifts the higher-frequency modes. 
To ensure the uos = 3uJi condition in the face of this 
dispersion, we slightly perturbed the structure (increas- 
ing the dielectric constant slightly at the nodes of the 
third-harmonic eigenfield) to tune the frequencies. The 
simulated crystal was effectively semi-infinite, with many 
more layers on the right than on the left of the cavity. On 
the left of the cavity, after two period of the crystal the 
material is simply air (e = 1), terminated by a perfectly 
matched layer (PML) absorbing boundary region. 

We excite the cavity with an incident plane wave of 
frequency cJi, and compute the resulting refiection spec- 
trum. The refiected power at c^s, the third-harmonic 
generation, was then compared with the prediction of 
the coupled-mode theory. The frequencies, decay rates, 
and a and P coefficients in the coupled-mode theory were 
computed from a linear FDTD simulation in which the 
eigenmodes were excited by narrow-band pulses. The 
freely available FDTD code of [57] was employed. 

The results are shown in Fig. |2j in which the out- 
put power at coi and cos = Scoi is denoted by 
and |s3_p, respectively, while the input power at uji 
is denoted by In particular, we plot conve- 

nient dimensionless quantities: the third-harmonic con- 
version efficiency |53_ |^ / as a function of the di- 
mensionless product 712 ki+l^ in terms of the standard 
Kerr coefficient n2 = Sx^^V^^^- There is clear agree- 
ment between the FDTD and GMT for small values of 
722 ki+l^ (in which limit the conversion goes quadrati- 
cally with 77,2 However, as the input power in- 
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FIG. 2: Log-log plot of |s3-|Vki+|^ vs. nslsi+l^ for 
the coupled- mode theory (grey) and FDTD (black squares), 
where 722 is being varied. Also shown are the corresponding 
An/72 (dashed blue) and Acji/cji (solid red) curves. 



creases, they eventually begin to disagree, marking the 
point where second-order corrections are required. This 
disagreement is not a practical concern, however, because 
the onset of second-order effects coincides with the lim- 
its of typical materials, which usually break down for 
An/n = x*^^^max|£;p/2£ > 1%. This is why we also plot 
the maximum index shift /S.n/n in the same figure. 

Also shown in Fig. [2] is a plot of Acji/cJi = Re[(5cJi/cc;i]. 
As expected, when Acji is of the order of l/Qi ~ 
10~^, the frequency shift begins to destroy the frequency 
matching condition, substantially degrading the third- 
harmonic conversion. (It might seem that /S.n/n and 
/S^uoi/uoi should be comparable, but this is not the case 
because /S.n/n is the maximum index shift while Acji/cji 
is due to an average index shift.) 

More specifically, the details of our simulation are as 
follows. To simulate a continuous wave (GW) source 
spectrum in FDTD, we employ a narrow-bandwidth 
gaussian pulse incident from the air region, which ap- 
proximates a GW source in the limit of narrow band- 
width. This pulse is carefully normalized so that the 
peak intensity is unity, to match the GMT. The field 
in the air region is Fourier transformed and subtracted 
from the incident field to yield the reflected flux. Using 
only two periods of quarter-wave stack on the left of the 
cavity we obtained two cavity modes with real frequen- 
cies ui = 0.31818 (27rc/a), UJ2 = 0.95454 (27rc/a) and 
quality factors Qi = 1286 and Qs = 3726, respectively. 
Given these fleld patterns, we computed the aij and Pi 
coeflicients. We obtained the following coupling coef- 
flcients, in units of x^^V^- ^11 — 4.7531 x 10~^,a22 = 
5.3306xl0-'^,ai2 = 0^21 = 2.7847x lO""^, A = (4.55985- 
0.7244z) X 10-^ 
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COMPLETE FREQUENCY CONVERSION 

We now consider the conditions under which one may 
achieve complete frequency conversion: 100% of the in- 
cident power converted to output at the second or third 
harmonic frequency. As we shah see, this is easiest to 
achieve in the x^^^ case, and requires additional design 
criteria in the x^^^ case. 

The key fact in a x^^^ medium is that there are no 
frequency-shifting terms {a = 0), so the resonance con- 
dition 002 = 2a;i is not spoiled as one increases the power. 
The only requirement that we must impose is that exter- 
nal losses such as absorption are negligible (re,/c ^ Ts^k)- 
In this case, 100% conversion corresponds to setting 
5i_ = in the steady-state. Using this fact, Eqs. |3]|4] 
for an input source = 5i+ exp(za;it) yields the fol- 

lowing condition on the input power for 100% conversion: 



u^lm^ Tsar? 



2|/3irQ2g? 



(19) 



(A similar dependence of efficiency on Q1Q2 was previ- 
ously observed [51, ^ , although a critical power was not 
identified.) Thus, we can always choose an input power 
to obtain 100% conversion. If Qi ~ Q2, then this criti- 
cal power scales as V/Q^ where V is the modal volume 
(recah that p - l/W). 

This is limited, however, by our first-order approxima- 
tion in perturbation theory: if the input power becomes 
so large that second-order effects (or material break- 
down) become significant, then this prediction of 100% 
conversion is no longer valid. The key condition is that 
the fractional change in the refractive index be small: 
An/n <C 1. This can always be satisfied, in principle: if 
one chooses Qi and/or Q2 to be sufficiently large, then 
the critical power can be made arbitrarily small in prin- 
ciple. Not only does the critical power decrease with Q^, 
but the field intensity in the cavity (|aip) decreases as 
V/QiQ2^ and thus one can avoid large An/n as well as 
lowering the power. (Note that the field intensity goes as 
l/Q^ while the power goes as l/Q^ simply because the 
energy and power are related by a time scale of Q.) 

To illustrate second-harmonic conversion for a x^^^ 
medium, we plot the solution to the coupled-mode equa- 
tions as a function of input power in Fig. d The 100% 
conversion at the predicted critical power is clearly vis- 
ible. For this calculation, we chose modal parameters 
similar to the ones from the FDTD computation before: 
uji 0.3, U2 = 0.6, Qi = 10"^, Q2 = 2 X 10"^, with dimen- 
sionless pi = (4.55985 - 0.7244) x 10"^ 

A x^^^ medium, on the other hand, does suffer from 
nonlinear frequency shifts. For example. Fig. |2j which 
is by no means the optimal geometry, exhibits a maxi- 
mal efficiency of |53_p/|5i+p ^ 4 x 10~^, almost three 
orders of magnitude away from complete frequency con- 
version. On the other hand, we can again achieve 100% 
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conversion if we can force aij = 0, which can be done 
in two ways. First, one could employ two x*^^^ materials 
with opposite-sign x*^^^ values (e.g., as in Ref. |58|). For 
example, if the x^^^ is an odd function around the cavity 
center, then the integrals for a^j will vanish while the /3 
integrals will not. (In practice, a <C /3 should suffice.) 
Second, one could pre-compensate for the nonlinear fre- 
quency shifts: design the cavity so that the shifted fre- 
quencies, at the critical power below, satisfy the resonant 
condition cos + Alu^ = 3(c<;i -|- Aooi). Equivalently, design 
the device for a^j = and then adjust the linear cavity 
frequencies a posteriori to compensate for the frequency 
shift at the critical power. (This is closely analogous to 
the cavity detuning used for optical bistability [55 , in 
which one operates off-resonance in the linear regime so 
that resonance occurs from the nonlinear shift.) 

If aij is thereby forced to be zero, and we can also 



neglect external losses (absorption, etc. 
100% third-harmonic conversion {si- - 
when: 



as above, then 
0) is obtained 



l^i+l 



^1 |A|^<1^^,3 



1/2 



m\ QiQs 



1/2 



(20) 

If Qi ^ Q3, then this critical power scales as V/Q^ where 
V is the modal volume (recall that /3 ~ ^/V). This is 
precisely the scaling that was predicted for the power to 
obtain nonlinear bistability in a single- mode cavity |59j. 
Similarly, one finds that the energy density in the cavity 
(|aip) decreases proportional to V j \/Q\Qz' 

We demonstrate the third-harmonic conversion for 
Oiij = by plotting the solution to the coupled-mode 
equations as a function of input power in Fig. [4j Again, 
100% conversion is only obtained at a single critical 
power. Here, we used the same parameters as in the 
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FDTD calculation, but with a = 0. In this case, compar- 
ing with Fig. |2| we observe that complete frequency con- 
version occurs at a power corresponding to An/n ^ 10~^. 
This is close to the maximum power before coupled- 
mode/perturbation theory becomes invalid (either be- 
cause of second-order effects or material breakdown), but 
we could easily decrease the critical power by increasing 
Q. 

For both the x*^^^ and the x*^^^ effects, in Figs, [sjlij we 
see that the harmonic conversion efficiency goes to zero 
if the input power (or x) is either too small or too large. 
It is not surprising that frequency conversion decreases 
for low powers, but the decrease in efficiency for high 
powers is less intuitive. It corresponds to a well-known 
phenomenon in coupled- mode systems: in order to get 
100% transmission from an input port to an output port, 
the coupling rates to the two ports must be matched 
in order to cancel the back-reflected wave [50l |60]. In 
the present case, the coupling rate to the input port is 
~ 1/Qi, and the coupling rate to the output "port" (the 
harmonic frequency) is determined by the strength of the 
nonlinear coupling. If the nonlinear coupling is either 
too small or too large, then the rates are not matched 
and the light is reflected instead of converted. (On the 
other hand, we find that for large input powers, while the 
conversion efficiency as a fraction of input power goes to 
zero, the absolute converted power (|s2-P or |53_p) goes 
to a constant.) 

Finally, let us consider one other potential problem. 
Any physical x^^^ medium will generally also have x^^^ 7^ 
0, so if the power is large enough this could conceiv- 
ably cause a frequency shift that would spoil the second- 
harmonic resonance even in the x^^^ device. Here, we 
perform a simple scaling analysis to determine when this 
will occur. (Although the frequency shifting could poten- 



tially be compensated for as described above, one prefers 
that it be negligible to begin with.) In order to pre- 
serve the resonance condition, any fractional frequency 
shift Aoo/oo must be much smaller than the bandwidth 
l/Q, or equivalently we must have QAco/co <C 1. From 
above, Au ~ a;a|ap, and |ap ~ \si-^\'^Q/u. Suppose 
that we are operating at the critical input power P*^^^ 
for second-harmonic conversion, from Eq. 19 It then 



follows that we desire QAujjuj ^ Q'^a/ujP^ < 1. It 
is convenient to re-express this relationship in terms of 
P^^) ijj I f5Q^ ^ the third-harmonic critical power from 
Eq. [19) by assuming a ~ as discussed in the previous 
section. We therefore find that x*^^^ self-phase modula- 
tion can be ignored for x^^^ second-harmonic generation 
as long as p(2)/p(3) <^ 1. As discussed in the conclud- 
ing remarks, this is indeed the case for common materi- 
als such as gallium arsenide, where p(2)/p(3) ~ 1/30 for 
Q ~ 1000 and for typical values of the cavity lifetime and 
volume. Moreover, since p^'^) j p^^) r^XjQ^ one can make 
the ratio arbitrarily smaller if necessary (at the expense 
of bandwidth) by increasing Q. 



THE EFFECT OF LOSSES 

In practice, a real device will have some additional 
losses, such as linear or nonlinear absorption and radia- 
tive scattering. Such losses will lower the peak conversion 
efficiency below 100%. As we show in this section, their 
quantitative effect depends on the ratio of the loss rate 
to the total loss rate l/Q. We also solve for the critical 
input power to achieve maximal conversion efficiency in 
the presence of losses. 

For a x^^^ medium with a linear loss rate l/re,/c, we 
solve Eqs |3||4] for 1^2- P and enforce the condition for 
maximal conversion efficiency: (j|gf^|2 (['^2- P/I^i+P) = 
0. We thus obtain the following optimal input power and 
conversion efficiency: 



2r, 



l^i+l 



1^2-1 



(21) 
(22) 



Eq. 



It immediately follows that for zero external losses, i.e. 

gives 100% conversion and Eq. 21 re- 
For small external losses Tg^k <^ ^e,ki 



duces to Eq. 19 



22 



the optimal efficiency is reduced by the ratio of the loss 
rates, to first order: 



1^2- 



l^l+l 



(23) 



(A similar transmission reduction occurs in coupled- 
mode theory when any sort of loss is introduced into a 
resonant coupling process [54^.) 
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The same analysis for x^^^ yields the following critical 
input power and optimal efficiency: 



|g3-| 

ki+l 



1/2 



T1T3 
Ts,lTs,3 



(24) 
(25) 



where by comparison with Eq. [22j a first-order expan- 
sion for low-loss yields an expression of the same form as 
Eq. [23j the efficiency is reduced by the ratio of the loss 
rates, with r2 replaced by ra. 

A x^^^ medium may also have a nonlinear "two- 
photon" absorption, corresponding to a complex-valued 
X*^^^ , which gives an absorption coefficient proportional to 
the field intensity. This enters the coupled-mode equa- 
tions as a small imaginary part added to a, even if we 
have set the real part of a to zero. (The corresponding 
effect on P is just a phase shift.) That yields a nonlinear 
(NL) Te^k of the following form, to lowest order in the 
loss: 



1 

IniT 

^e,l 



^e,3 



: UJi Im 



! LOs Im 



2 3 

'^s,! I |2 I '^s,3'^s,l 2 1 /o |2 1 16 



(26) 



2 3 

I |2 I ^s,3^s,l 2\n \2\ 16 

«31^|S1+| +^33 ' ^z\H kl+l 



(27) 



where we have simply substituted the values for the crit 
ical fields a\ 



Y^2/ri5i+ and as given by Eq. [6j and 
grouped terms that correspond to imaginary frequency 
shifts. These loss rates can then be substituted in the 
expression for the losses above, i.e. Eq. [25) in which 
case one obtains the following optimal efficiency of third- 
harmonic generation, to lowest-order in the loss, not in- 
cluding linear losses: 



1^3- 



l^l+l 




ail + 3ai3 , ai3 + 3^33 



Ts,3 



(28) 

(The linear and nonlinear losses can be combined by sim- 
ply multiplying Eq. 25 and Eq. 28 ) Thus, the nonlinear 
loss is proportional to the ratio Ima/|/5|, which is pro- 
portional to Imx'^^Vlx'^^^l- 



CONCLUSION 

We have presented a rigorous coupled-mode theory for 
second- and third-harmonic generation in doubly reso- 
nant nonlinear cavities, accurate to first order in the non- 
linear susceptibility and validated against a direct FDTD 
simulation. Our theory, which generalizes previous work 



on this subject, predicts several interesting consequences. 
First, it is possible to design the cavity to yield 100% 
frequency conversion in a passive (gain- free) device, even 
when nonlinear down- conversion processes are included, 
limited only by fabrication imperfections and losses. Sec- 
ond, this 100% conversion requires a certain critical input 
power — powers either too large or too small lead to lower 
efficiency. Third, we describe how to compensate for the 
self-phase modulation in a x^^^ cavity. The motivation 
for this work was the hope that a doubly resonant cavity 
would lead to 100% conversion at very low input powers, 
and so we conclude our paper by estimating the critical 
power for reasonable material and geometry assumptions. 



A typical nonlinear material is gallium arsenide 
(GaAs), with x^^^ - 145 pm/V and = 1.5 x 10"^^ 
cm^/W at 1.5/im. (Al doping is usually employed to de- 
crease nonlinear losses near resonance [61 .) Although 
this has both x*^^^ and x*^^^ effects, we can selectively en- 
hance one or the other by choosing the cavity to have 
resonances at either the second or third harmonic. Many 
well confined optical cavity geometries are available at 
these wavelengths and have been used for nonlinear de- 
vices, such as ring resonators [62 or photonic-crystal 
slabs [63]. We will assume conservative parameters for 
the cavity: a lifetime Qi = 1000, Q2 = 2000, Q3 = 3000, 
and a modal volume of 10 cubic half-wavelengths {V ~ 
10(A/2n)^) with roughly constant field amplitude in the 
nonlinear material (worse than a realistic case of strongly 
peaked fields). In this case, the critical input power, from 
Eqs. [l9]-[20) becomes approximately 70 /iW for second- 
harmonic generation and 2 mW for third-harmonic gen- 
eration (with a moderate peak index shift An/n 10~^, 
justifying our first-order approximation) 



Future work will involve designing specific doubly res- 
onant cavity geometries and more precise power predic- 
tions. Using our expressions for a and /3, optimized cav- 
ities for harmonic generation can be designed using stan- 
dard methods to compute the linear eigenmodes. In prac- 
tice, experimentally achieving cavity modes with "ex- 
actly" harmonic frequencies, matched to within the frac- 
tional bandwidth l/Q, is a challenge and may require 
some external tuning mechanism. For example, one could 
use the nonlinearity itself for tuning, via external illumi- 
nation of the cavity with an intense "tuning" beam at 
some other frequency. Also, although we can directly in- 
tegrate the coupled-mode equations in time, we intend 
to supplement this with a linearized stability analysis at 
the critical power. This is particularly important for the 
X^^^ case, where pre-correcting the frequency to compen- 
sate the nonlinear frequency shift (self-phase modulation) 
may require some care to ensure a stable solution. 
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